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Executive  Summary 

This  report  covers  the  activities  at  the  University  of  Notre  Dame  under  AFOSR  Grant 
F49620-97- 1-0489,  “Fluid-Optic  Interactions  II.”  The  report  begins  with  a  discussion 
of  the  relevance  of  the  research  to  Air  Force  missions  and  places  the  work  in  context 
with  work  performed  prior  to  the  time  the  grant  effort  was  begun.  When  a  laser 
beam  propagates  through  a  variable-index-of-refraction,  turbulent  fluid,  its  optical 
wavefront  becomes  aberrated,  reducing  associated  optical  system  performance.  At 
flight  Mach  numbers  above  ~  0.6,  “compressibility  effects”  alone  become  important 
in  aberrating  wavefronts.  Scaling  optical  wavefront  data  collected  in  ground-test  facil¬ 
ities  to  realistic  flight  conditions  requires  an  understanding  of  how  this  compressibility 
mechanism  produces  wavefront  distortions.  Prior  to  this  grant’s  research,  little  was 
known  about  this  compressibility  mechanism.  This  compressibility  mechanism  was 
studied  using  an  analytical/ numerical  model.  Because  the  shear-layer  flows  in  ques¬ 
tion  have  low  convective  Mach  numbers,  they  can  be  considered  weakly-compressible, 
which  allows  their  velocity  fields  to  be  approximated  using  a  discrete  vortex  method. 
From  these  “known”  velocity  fields,  the  concomitant  pressure  and  density  fields  were 
determined  by  iteratively  solving  the  unsteady  Euler  equations.  The  results  showed 
that  unsteady  pressure  fluctuations  present  in  shear  layers  had  a  dominant  effect  on 
the  resulting  density /index  fields;  this  effect  has  historically  been  neglected.  The 
computed  index  fields  produced  simulated  schlieren  images  which  closely  resembled 
experimental  schlierens.  Optical  wavefronts  computed  from  the  simulation  reasonably 
matched  the  behavior  of  large-scale  aberrations  measured  in  a  transonic  wind  tunnel 
at  AEDC.  An  explanation  for  small-scale  distortions  observed  in  these  experimental 
data  was  formulated,  suggesting  that  temperature  discontinuities  in  the  splitter-plate 
boundary  layer  had  been  fed  into  the  shear  layer.  The  compressibility  model  was  used 
to  suggest  and  test  relations  for  scaling  optical  wavefronts  from  one  flow  condition  to 
another.  A  simple  density  ratio  successfully  scaled  distortion  variations  with  altitude 
when  shear  layer  Mach  numbers  were  held  constant;  an  additional  temperature  ratio 
was  required  if  the  velocity  difference  was  kept  the  same.  Optical  distortions  produced 
by  dissimilar-index,  incompressible,  mixing  flows  were  found  to  have  a  different  char¬ 
acter  than  those  from  weakly-compressible  flows;  however,  the  spatial  and  temporal 
frequencies  of  low-speed  distortions  were  found  to  logically  scale,  provided  the  beam’s 
diameter  was  larger  than  the  largest  flow  structures.  The  results  of  the  research  were 
communicated  directly  to  DoD  personnel  working  in  the  technical  area  as  well  as 
being  presented  at  several  national  technical  meetings.  The  work  is  documented  in  a 
Ph.D.  Dissertation,  six  meeting  papers,  and  four  journal  articles. 
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Final  Report 

FLUID-OPTIC  INTERACTIONS  II 
AFOSR  Grant  F49620-97- 1-0489 

Principal  Investigator:  Eric  J.  Jumper 
Hessert  Center  for  Aerospace  Research 
University  of  Notre  Dame 
Notre  Dame,  Indiana,  46556,  USA 
(219)  631-7680 
FAX:  (219)  631-8355 
Email:  Eric.J.Jumper.l@nd.edu 

1.  INTRODUCTION 

The  use  of  optical  signals  propagated  through  fluids  is  ubiquitous  in  Air  Force  and 
DoD  applications.  These  applications  include,  but  are  not  limited  to,  imaging  distant 
objects  from  ground-based  and  airborne  systems,  tracking,  communication,  ranging, 
sensing  incoming  threats,  irradiating  targets  for  weapons-delivery  identification,  and 
destroying  targets  by  irradiation  with  high-power  lasers.  These  latter  systems  include 
airborne  laser  weapons  such  as  the  Air  Force’s  Airborne  Laser  (ABL)  theater  ballistic 
defense  system. 

Aberrations  (distortions  of  figure)  on  optical  wavefronts  degrade  the  ability  of 
an  optical  system  to  image,  focus,  or  otherwise  use  the  optical  signal  associated  with 
the  propagating  wavefront  for  the  system’s  intended  purpose.  When  an  otherwise- 
planar  (flat-figured)  wavefront  propagates  through  an  optically-active  field  (i.e. ,  a 
field  of  variable  index  of  refraction),  the  wavefront  emerges  aberrated.  If  the  im¬ 
posed  aberrations  are  known  and  constant,  measures  may  be  taken  to  compensate 
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for  the  aberrations  to  improve  the  performance  of  the  optical  system;  however,  when 
the  aberrating  field  is  rapidly  time-varying,  the  problem  of  compensation  becomes 
more  complicated.  In  the  case  of  a  fluid-optic  interaction,  the  distortion  is  caused  by 
a  dynamically-changing  turbulent  flowfield;  thus,  the  wavefront’s  aberrations  them¬ 
selves  are  also  dynamically  changing.  Although  the  spatial  and  temporal  frequencies 
of  the  optical  aberrations  are  related  to  the  spatial  and  temporal  character  of  the  tur¬ 
bulent  flowfield,  the  optical  aberrations  result  from  an  integrated  path  through  the 
variant-index  field  associated  with,  but  not  identical  to,  the  turbulent  flowfield.  As 
such,  optical  aberrations  cannot  be  directly  related  to  any  particular  instantaneous 
spatial  feature  in  the  flowfield;  however,  the  fact  that  the  aberrations  are  associated 
with  features  in  the  turbulent  flowfield,  features  which  convect  at  some  velocity  in 
the  flow  direction,  means  that  the  aberrated  features  on  the  wavefront  also  “convect” 
in  the  flow  direction.  This  fact  has  been  of  critical  importance  to  our  research. 

1.1.  Defining  the  Fluid-Optic  Interaction  Problem 

This  system-performance  reduction  induced  by  turbulent,  variant-index  flowfields 
can  be  quantified  by  analysis  of  the  far-field  irradiance  pattern.1  Light  transmis¬ 
sion  through  such  time-varying,  aberrating  fields  can  occur  over  short  path  lengths 
in  the  “near  field”  (i.e.,  path  length,  y,  on  the  order  of  the  system  aperture  diameter, 
A),  or  can  occur  over  extended  paths  (i.e.,  path  lengths,  y,  much  larger  than  A).  It 
should  be  noted  that  the  term  “near  field”  in  the  present  context  refers  to  aberrating 
fluid  structures  very  near  the  exit  aperture  of  an  outgoing  collimated  beam,  or  at 
the  receiving  aperture  for  an  imaging  system.  This  should  not  be  confused  with  the 
“near  field”  in  physical  optics,  where  the  Fresnel  approximations,  but  not  Fraunhofer 
approximations  can  be  applied.  The  former,  short-path  (near-field)  aberration  is 
usually  termed  “aero-optics”  and  the  latter,  long-path  aberrations  are  usually  termed 
“atmospheric  propagation.”  These  classifications  have  their  origin  in  the  specific  ap- 
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plications  with  which  they  are  associated  (see  Refs.  2  and  3,  for  example).  It  has 
long  been  known  that  placing  a  conjugate  waveform  on  the  beam  prior  to  its  trans¬ 
mission  through  the  aberrating  medium  (or  prior  to  focusing  for  a  system  receiving 
an  optical  signal  that  has  traversed  the  aberrating  medium)  results  in  the  emergence 
of  a  planar-wavefront  beam  as  it  leaves  the  medium  (or  enters  a  focusing  mechanism 
after  traversing  the  medium);  systems  that  sense  the  aberration,  construct  and  apply 
the  proper  conjugate  waveform  at  regular  time  intervals  are  termed  “adaptive-optic 
systems,”  and  the  study  of  such  systems  is  known  as  “adaptive  optics.”4  The  deter¬ 
mination  of  the  design  requirements  for  such  systems  depend  on  understanding  the 
spatial  and  temporal  frequencies  of  the  distorted  optical  wavefront  for  the  applicable 
aberrating  flowfield. 

Over  the  last  decade,  progress  has  been  made  in  both  measuring  the  wave- 
front  dynamics  for  the  atmospheric  propagation  problem  and  using  this  information 
for  designing  and  applying  adaptive-optic  equipment  and  techniques.5  For  aero-optic 
problems,  on  the  other  hand,  progress  in  adaptive-optic  correction  has  been  essentially 
nonexistent.  The  reason  for  this  is  that  the  required  spatial  and  temporal  frequencies 
associated  with  aero-optic  problems  are  at  least  an  order  of  magnitude  greater  than 
those  presently  correctable  by  adaptive-optic  systems  for  the  atmospheric-propagation 
case,  even  for  relatively-slow  laboratory  flows.5,6  As  such,  the  main  research  thrust  in 
aero-optics  until  very  recently  has  been  to  attempt  to  quantify  the  time- averaged,  sta¬ 
tistical,  near-field  wavefront  distortion  (as  phase  variance  or  root-mean-square  Optical 
Path  Difference,  OPDrms )  in  order  to  estimate  the  time-averaged,  far-field  degrada¬ 
tion  that  might  be  expected  for  a  system  whose  optical  signals  must  traverse  the 
aberrating,  turbulent  flowfield.7  In  this  context,  “far  field”  is  defined  as  the  region 
where  the  diffraction  integral  can  be  determined  using  Fourier  Optics. 

It  is  helpful  in  understanding  fluid-optic  interaction  problems  to  divide  the 
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problems  into  four  levels  of  understanding  associated  with  either  the  atmospheric 
propagation  problem  or  the  aero-optic  problem.  These  levels  may  be  generally  ar¬ 
ranged  in  their  order  of  complexity  as  follows:  first,  there  is  the  level  of  understanding 
that  allows  one  to  estimate  the  statistical  level  of  optical  distortion  that  is  likely  to  be 
encountered  for  a  particular  optical  propagation  scenario  through  a  flowfield.  At  the 
second  level,  one  can  predict  not  only  the  statistically-averaged  distortion  but  also 
the  spatial  and  temporal  frequencies  associated  with  the  time-varying  optical  distor¬ 
tion.  The  third  level  is  reached  when  actual  time  histories  associated  with  the  optical 
distortions  can  be  measured  and  linked  to  the  specific  fluid-mechanic  character  of  the 
flowfield.  With  a  fourth-level  understanding,  the  optical  distortions  can  be  measured 
rapidly  enough  to  allow  adaptive-optic  compensatation  for  the  distortions,  thereby 
improving/restoring  the  optimal  performance  of  the  optical  system.  The  atmospheric 
propagation  problem  has  advanced  through  all  four  levels  of  understanding,  while 
until  the  span  of  time  that  covers  our  initial  AFOSR  funding,  the  aero-optic  problem 
remained  at  level  one,  and  the  researchers  in  the  field  had  no  prospects  or  apparent 
aspirations  to  move  to  a  higher-level  understanding.  A  historical  review  of  aero-optics 
is  given  as  APPENDIX  A  that  attempts  to  place  the  statistical  approach  in  perspec¬ 
tive.  The  review  includes  reference  to  Refs.  8-9  that  span  the  period  from  19528  to 
1992. 9  Our  first  grant  (F49620-93- 1-0163)  began  one  year  later  in  1993. 

1.2.  Association  of  Aero-Optics  with  Fluid  Mechanics. 

It  is  clear  from  APPENDIX  A  that  the  aero-optic  problem  has  always  been  inextri¬ 
cably  linked  with  the  fluid  mechanics  of  turbulent  flows.  This  is  in  contrast  to  the 
atmospheric-propagation  problem  which  has  been  treated  almost  exclusively  as  an 
optical  problem.  For  astronomical  observations,  assuming  that  viewing  conditions 
are  favorable,  the  aberrating  structures  are  large  compared  to  the  viewing  aperture. 
Often  these  structures  are  of  such  extent  that  the  aberrations  can  be  treated  as  sim- 
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pie  tilt  corrections,5  and  smaller  aberrations  are  still  relatively  large  with  respect 
to  the  aperture,  so  that  only  the  first  few  aberration  modes  need  be  addressed  in 
order  to  restore  much  of  the  aberration-degraded  system  performance.  These  facts 
mean  that  the  resolution  of  a  wavefront  sensor  to  detect  the  aberrations  need  not 
be  exceptionally  high.  Beyond  this,  the  aberrations  “convect”  through  the  viewing 
aperture  at  speeds  associated  with  atmospheric  winds  plus  any  slue  rates  that  may  be 
superimposed.  Typical  rates  required  to  provide  time-resolved  measurements  of  the 
aberrated  wavefronts  due  to  propagation  through  the  atmosphere  are  on  the  order  of 
hundreds  of  Hertz.  Wavefront  sensor  technology  has  for  some  time  been  capable  of 
achieving  these  rates.4,5 

The  aero-optic  problem,  by  its  very  nature,  must  contend  with  shear-layer  and 
boundary-layer  aberrations  whose  sizes  are  on  the  order  of  the  associated  turbulent 
structures;  further,  the  convection  rates  are  now  associated  with  the  defining  flows 
which  can  be  much  higher  than  typical  atmospheric  winds  (in  the  case  of  flow  over  an 
aircraft  in  flight,  for  example).  Our  efforts  have  shown  that  even  for  a  laboratory  flow 
of  less  than  10  m/s,  because  of  the  small  size  of  the  turbulent  structures,  wavefront 
capture  rates  of  5  kHz  are  needed  to  time-resolve  the  dynamics  of  the  aero-optic 
wavefront  aberrations.  Until  very  recently,  the  fastest  wavefront  sensors  available 
operated  at  less  than  1  kHz  (see  discussion  of  Shack-Hartmann  Sensors  in  APPENDIX 
B).  For  this  reason,  equations  linking  the  statistical,  optical-aberration  field  over 
the  aperture,  like  Eq.  12  in  APPENDIX  A,  were  developed.  Examination  of  Eq.  12 
shows  that  this  estimate  depends  entirely  on  statistical  measures  of  the  fluid-mechanic 
properties  of  the  turbulent  flowfield  over  the  viewing  aperture. 

1.3.  The  Aero-Optic,  Adaptive-Optic  Problem 

As  mentioned  above,  in  order  to  even  scope  the  order  of  complexity  associated  with 
attempting  to  provide  an  adaptive-optic  correction  for  optical  propagation  through 
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an  aero-optic  field,  a  level-two  understanding  of  the  aberrations  is  needed.  In  1992, 
practitioners  in  the  field  of  aero-optics  felt  that,  with  no  hope  of  a  significant  break¬ 
through  in  wavefront  sensor  technology,  the  field  of  aero-optics  could  never  mature 
to  the  level  of  addressing  adaptive-optics;  however,  some  few  of  us  were  at  least  curi¬ 
ous  about  the  requirements  of  an  adaptive  optic  system  for  the  aero-optics  problem. 
My  early  attempts  at  determining  estimates  for  the  spatial  and  temporal  frequencies 
relevant  to  the  aero-optics  problem  using  the  statistical  measures  of  turbulence  were 
less  than  satisfactory,  and  served  as  the  basis  for  our  first  AFOSR  grant.  That  grant 
initially  sought  only  to  address  this  question;  i.e.,  it  sought  to  obtain  a  level  two 
understanding  of  the  aero-optic  problem.  Work  performed  in  the  late  1980’s,  had 
suggested  to  us  that  there  might  be  a  way  of  measuring  these  spatial  and  temporal 
frequencies.  In  Refs.  10  and  11,  Malley  et  al.  describe  a  simple  optical  device  that, 
when  coupled  with  certain  fluid-mechanic  principles,  was  capable  of  directly  measur¬ 
ing  the  statistical  OPDrms.  APPENDIX  B  describes  the  principles  associated  with 
the  instrument,  along  with  some  background  on  Hartmann  wavefront  sensors.12 

In  that  first  round  of  funding,  we  soon  came  to  realize  that  the  fundamental 
assertion  upon  which  the  Malley  et  al.  instrument  was  based  had  the  potential 
of  developing  an  actual  wavefront  sensor.  One  major  accomplishment  of  the  first 
grant  was  the  development  of  the  SABT  (Small-Aperture-Beam  Technique)  wavefront 
sensor.6,13  The  (one-dimensional)  SABT  sensor  is  described  in  APPENDIX  C.  In  this 
round  of  funding,  a  two-dimensional  version  of  the  sensor  was  developed  as  described 
below.14,15  The  SABT  sensor  made  possible  for  the  first  time  the  ability  to  examine 
the  time-resolved  dynamics  of  the  aero-optically  aberrated  wavefront.  We  studied 
wavefronts  distorted  by  propagation  through  turbulent  shear  flows  between  fluids 
of  dissimilar  index  of  refraction;  in  particular,  we  extensively  studied  propagation 
through  the  transitional  region  of  a  heated,  two-dimensional  jet.  Fig.  1,  produced 
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using  a  numerically-simulated  flowfield,1"^  is  helpful  in  understanding  the  nature  of  the 
optical  aberrations  produced  by  propagation  through  the  heated  jet.  At  an  instant 
in  time,  an  initially-planar  wavefront  propagates  through  the  aberrating  flowfield, 
also  shown  at  that  instant,  and  emerges  with  the  distorted  wavefront  shown.  In 
this  case,  the  aberration  is  due  to  the  mixing  of  dissimilar  index-of-refraction  flows. 
Although  the  distortion  is  due  to  a  path  integration  through  the  flow,  it  is  possible 
to  associate  features  on  the  wavefront  with  flow  features.  As  the  flowfield  convects 
and  evolves,  a  footprint  of  the  character  of  the  evolving  flowfield  is  captured  by  the 
wavefront;  a  time  series  of  the  instantaneous  wavefronts,  as  shown  at  the  upper  portion 
of  Fig.  1,  captures  the  dynamics  of  the  shear  layers  of  the  two-dimensional  jet.  The 
characteristic  hills  and  valleys,  and  the  slopes  of  these  features  with  time,  indicate 
the  convective  and  evolving  nature  of  the  flowfield;  different  structures  convect  at 
different  rates  as  evidenced  by  the  different  slopes  of  the  hills  and  valleys.  The  fact 
that  the  shear-layer  roll  up  is  organized  can  be  seen  by  the  well-ordered  nature  of  the 
time  series  of  wavefronts. 

An  actual  time  series  of  experimentally  measured  wavefronts  from  the  heated- 
jet  experiments  is  shown  in  Fig.  2.6  It  is  clear  that  the  features  found  in  the  simulated 
field  are  present  in  the  experiments.  Time  series  of  the  Fig.  2  type  were  studied  in 
order  to  assess  the  adaptive-optic  correction  requirements  that  would  be  needed  to  im¬ 
prove/restore  the  Strehl  ratio,  a  measure  of  performance  (see  APPENDIX  A,  Eq.  13). 
Fig.  3  gives  a  collage  of  representative  results  from  Ref.  16;  the  equation  defining  the 
Strehl  Ratio,  SR,  is  given  in  the  lower  left  corner  of  the  figure.  The  upper-right  corner 
of  Fig.  3  shows  the  far-field  focus  for  a  periodic  application  of  a  perfect  correction 
restoring  the  Strehl  ratio  to  1.0  each  time  it  is  applied.  In  between  these  applications, 
the  focus  degrades  as  the  aberration  changes  with  the  convecting  and  evolving  tur¬ 
bulent  structures.  In  this  particular  example,  the  uncorrected,  time-averaged  Strehl 
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Figure  1.  Schematic  of  a  planar  wavefront  propagating  through  a  two-dimensional 
heated  jet  and  concomitant  time  series  of  one-dimensional,  aberrated  wavefronts. 

ratio  was  0.265;  i.e.,  the  peak  intensity  of  the  focal  point  was  approximately  26%  of  a 
diffraction-limited  peak  intensity.  For  this  example,  propagation  through  a  7  m/s  jet, 
a  perfect  correction  would  have  to  be  applied  at  10  kHz  to  restore  the  Strehl  to  greater 
than  95%  diffraction-limited  intensity.  With  a  four-by-four-element  deformable  mir¬ 
ror  over  the  aperture  shown,  also  cued  at  10  kHz,  the  Strehl  ratio  would  be  restored 
to  approximately  75%  of  diffraction-limited.  The  lower-right  figure  indicates  that 
the  restoration  will  be  further  reduced  if  the  correction  applied  is  itself  “late.”  This 
means  that  in  order  to  maximize  the  restored  intensity  the  applied  correction  must 
be  correct  at  the  time  it  is  applied. 

Most  future  aero-optic  research  will  be  done  in  ground-test  facilities.  The  dif¬ 
ficulties  of  testing  in  compressible-flow  facilities  (cost,  test-section  size,  test  duration, 
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Figure  2.  Space-time  representation  of  SABT-measured  OPDs  for  propagation 


through  an  experimental,  two-dimensional,  heated  jet  flowfield  (from  Ref.  6). 


ADAPTIVE-OPTIC  CORRECTION 
WITH  A  FINITE  UPDATE  RATE 


FAR-FIELD  FOCUS 


SR  =  f(x1)g(x2) 

x,  =  time  between  updates 

x2  =  time  needed  to  sense 
and  construct 
wavefront  correction 


Figure  3.  Adaptive-optic  correction  requirements  study. 
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vibration)  will  limit  the  number  and  success  of  such  tests,  in  many  instances,  aero- 
optical  distortions  have  been  and  will  continue  to  be  produced  by  low-speed  (~  7  m/s, 
say)  mixing  of  dissimilar-index,  incompressible  flows.  The  presumption  inherent  in 
using  low-speed  flows  for  aero-optic  research  is  that  the  results  of  such  research  will 
also  be  applicable  to  high-speed  aero-optics.  The  question,  then,  is  how  would  wave- 
front  data  collected  in  ground-test  facilities  scale  to  other  flight  conditions  and,  in 
particular,  can  data  obtained  in  low-speed  flows  be  used  to  infer  distortions  created 
by  high-speed  flows? 

In  the  last  stages  of  our  first  funding  round,  we  used  the  SABT  sensor  to 
measure  the  wavefront  dynamics  of  propagation  through  the  Mach  0.8  compressible 
shear  layer  at  the  Arnold  Engineering  Development  Center  (AEDC). 17,18  Fig.  4  shows 
an  example  time  series  of  experimental  wavefronts  from  these  tests;  these  data  were 
taken  at  station  2,  roughly  one  half  meter  downstream  from  the  splitter  plate.  It 
is  clear  from  the  figure  that  the  dynamics  of  the  wavefront  aberrations  are  indeed 
similar  in  character  to  those  for  propagation  through  the  heated  jet  (compare  Fig.  4 
to  Fig.  2,  for  example).  For  the  wavefront  structures  shown  in  Fig.  4,  when  compared 
to  those  in  Fig.  2,  the  sizes  of  the  AEDC  structures  (i.e.,  spatial  scale)  are  roughly 
eight  times  larger  than  those  wavefront  structures  for  the  heated  jet.  At  the  same 
time,  the  convection  velocity  of  the  structures  in  the  AEDC  experiments  increased 
by  roughly  forty  times.  Thus,  one  might  conclude  that  for  an  adaptive-optic  system 
to  perform  at  the  same  level  of  restoration  as  for  the  heated  jet,  using  a  similar 
number  of  subapertures,  would  require  corrections  to  be  made  at  60  kHz;  however, 
our  initial  attempts  at  scaling  heated-jet  OPD  structure  sizes/frequencies  to  predict 
the  AEDC  measurements  was  unsuccessful,  casting  doubt  on  how  adaptive-optics 
requirements  would  scale  for  other  system  apertures  and  conditions.  This  is  not  just 
an  academic  exercise.  A  large-scale  weapon  system  with  ~  1-m  aperture  is  being 
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prepared  for  fielding  (i.e.,  the  ABL);  yet,  there  was  no  method  for  estimating  aero- 
optical  effects  the  system  might  expect  prior  to  the  presentation  of  some  of  this 
grant’s  research  at  the  Second  Annual  Directed  Energy  Symposium  (1-5  November 
1999). 19,20  Thus,  scaling  issues  were  the  focus  of  a  large  portion  of  the  work  of  this 
round  of  funding  (F49620-97-0489).  The  individual  objectives  of  the  present  grant’s 
research  are  discussed  in  the  next  section. 


Figure  4.  AEDC  station  2  SABT-measured  wavefronts:  beam  separation  =  2.5  cm, 
High-Pass  Filter  =  2500  Hz.  (From  Ref.  18.) 

2.  OBJECTIVES  OF  THE  RESEARCH 

The  primary  objective  of  the  present  grant’s  research  was  to  address  the  scaling  issues 
that  came  to  light  in  studying  the  similarities  between  optical  aberrations  due  to  prop¬ 
agation  through  the  heated  jet  and  propagation  through  the  compressible  shear  layer. 
Additional  objectives  were  to:  (a)  extend  some  of  the  instrumentation-development 
work  of  the  last  grant  by  developing  a  derivative,  two-dimensional  wavefront  sensor  of 
the  one-dimensional  Small- Aperture  Beam  Technique  (SABT)  sensor;  (b)  investigate 
the  possible  use  of  control  to  develop  a  scheme  for  real-time  sensing  and  correcting  the 
wavefront  for  propagation  through  compressible  shear  layers;  and  (c)  further  explore 
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the  relationship  between  the  wavefront  aberrations  and  the  concomitant  aberrating 
structures.  The  approach  and  results  for  the  scaling  objective  are  described  in  Sec¬ 
tion  3.  The  approaches  and  results  for  the  remaining  objectives  are  described  in 
Sections  4-6. 

3.  SCALING  OF  OPTICAL  ABERRATIONS 
3.1.  Approach 

To  develop  reliable  scaling  relations  requires  a  basic  understanding  of  the  phenomenon 
to  be  scaled  (in  this  case,  compressibility-caused,  optical  aberrations).  Prior  to  this 
grant,  we  tried  scaling  the  spatial  and  temporal  frequencies  of  distortions  created  by  a 
low-speed,  heated  jet  using  relations  from  the  fluid-mechanics  literature,  but  we  could 
not  match  the  frequencies  measured  in  a  compressible  shear  layer  using  the  SABT 
sensor.21  Without  a  fundamental  understanding  of  the  compressibility  mechanism, 
the  reasons  for  this  difficulty  were  not  understood. 

The  mechanism  by  which  light  is  aberrated  by  transmission  through  a  shear 
layer  composed  of  two,  dissimilar-index  fluids  is  fairly  well  understood.22  Prior  to  the 
present  grant’s  research,  the  same  could  not  be  said  for  a  shear  layer  where  the  index- 
of-refraction  variations  are  due  to  compressibility.  Cassady’s23  Navier-Stokes  simula¬ 
tion,  for  instance,  only  modeled  aero-optical  distortions  produced  by  the  dynamics  of 
shear-layer  reattachment,  not  those  caused  by  shear  layer  coherent  structures.  Direct 
numerical  simulations  of  compressible  shear  layers  had  been  performed  by  the  early 
1990’s  24,25  but  these  studies  were  primarily  interested  in  shear-layer  stability  and  the 
resulting  coherent  structure  shapes  rather  than  on  the  instantaneous  density  fields  re¬ 
quired  for  aero-optical  calculations.  As  late  as  1998,  some  researchers  were  implicitly 
assuming  the  compressibility  mechanism  to  be  the  same  as  that  for  dissimilar-index 
mixing;  for  example,  Dimotakis  et  al.26  measured  optical  distortions  produced  by 
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a  compressible  mixing  layer  using  streams  of  markedly  different  index-of-refraction, 
thereby  overwhelming  any  aberrations  created  by  compressibility  alone.  The  lack  of 
compressibility-mechanism  understanding  became  apparent  to  us  when  the  first  com¬ 
pressibility  model  attempted  during  the  present  research  was  explained  to  a  group  of 
experts  in  compressible  shear  layers  and  in  aero-optics  at  a  1998  technical  meeting. 
While  all  in  attendance  agreed  that  this  initial  model  produced  unrealistically-low 
aberrations,  no  one  was  sure  what  was  missing  from  the  model.27-29 

Thus,  the  primary  goal  of  the  research  under  this  grant  was  to  delve  into  this 
question  of  the  compressibility  mechanism  and,  based  on  those  results,  suggest  how 
ground  (or  flight)  test  data  might  be  scaled.  As  such,  this  research,  in  effect,  consti¬ 
tutes  an  enabling  technology  required  for  future  success  in  adaptive-optic  correction 
of  aero-optic  distortions. 

The  present  research  began  with  two  things  that  could  be  exploited  from  our 
first  grant  (F49620-93-0163).  The  first  of  these  was  the  database  of  optical  wavefront 
measurements  serendipitously  obtained  in  a  compressible  shear  layer.  Although  some 
of  this  data  had  already  been  reported  in  the  literature,17,18  much  of  the  database 
had  not  been  examined  prior  to  the  present  grant.  This  wavefront  data  was  used  to 
anchor  the  analytical/numerical  model  developed  during  the  present  research. 

The  second  item  that  was  used  was  Hugo’s  numerical  shear-layer  simulation.22 
Hugo  developed  the  simulation  as  a  tool  for  testing  algorithms  required  in  the  devel¬ 
opment  of  the  SABT  wavefront  sensor.  For  his  purposes,  the  simulation  only  needed 
to  provide  a  deterministic,  index-of-refraction  field  that  evolved  as  it  convected  down¬ 
stream.22  The  present  study,  however,  further  required  that  the  simulation  compute  a 
realistic,  first-order  estimate  of  the  time- varying  velocity  field  produced  by  a  weakly- 
compressible  shear  layer.  These  velocity  fields  provided  a  common  base  flow  for 
further  analysis  and  allowed  direct  comparison  of  the  index  fields  and  resulting  aero- 
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optic  distortions  produced  by  both  incompressible  and  weakly-compressible  flows. 

Analytical/numerical  methods  were  then  used  to  determine  index-of-refraction 
(or  equivalently,  density)  fluctuations  in  a  weakly-compressible  shear  flow.  A  com¬ 
pressibility  model  was  developed  in  which  a  series  of  mechanisms  were  overlaid  on  the 
common  velocity  field  to  compute  the  time- varying  density  field.  Optical  wavefronts 
were  then  computed  by  (numerically)  propagating  an  optical  beam  through  the  time- 
varying  index  (density)  field.  The  validity  of  the  overlaid  mechanisms  was  assessed  by 
comparing  the  character  of  the  computed  aero-optical  distortions  to  those  measured 
experimentally  in  the  similar,  weakly-compressible,  shear  layer  discussed  above.18  The 
experimentally-anchored  compressibility  model  was  then  accepted  as  the  basis  for  de¬ 
veloping  appropriate  relations  for  scaling  aero-optic  distortions  to  other  conditions 
(e.g.,  ground  test  to  flight). 

3.2.  Re-reduction  of  AEDC  Experimental  Data 

The  AEDC  optical  wavefront  data,  first  presented  and  examined  in  Refs.  17  and  18, 
documented  the  optical  distortions  created  by  a  compressible  shear  layer  over  a  lim¬ 
ited,  5-cm,  test  aperture.  The  AEDC  test  section  is  shown  schematically  in  Fig.  5. 
As  discussed  in  Ref.  17,  a  post- test,  high-pass  digital  filter  was  employed  to  remove 
signal  noise  caused  by  vibration  of  the  test  section’s  optical  access  windows.  In  this 
earlier  work,  the  cutoff  frequency  was  conservatively  chosen  in  order  to  minimize  the 
chances  of  vibrational  corruption  of  the  beam  jitter  data;17  the  2-2.5-kHz  corner  fre¬ 
quency  was  also  chosen  because  the  small  measurement  aperture  already  functioned 
as  a  spatial  filter  for  any  larger-scale  structures.  For  this  data  to  be  used  to  anchor 
the  analytical/numerical  models  developed  during  this  grant,  a  lower  corner  frequency 
was  desired  in  order  to  retain  the  influences  of  flow  structures  larger  than  the  test 
aperture,  as  explained  below  and  documented  in  Refs.  30-32. 
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Figure  5.  Schematic  of  AEDC  facility  (From  Ref.  33). 

The  growth  rate  of  a  shear  layer  is  now  well  understood  to  be  driven  by  the 
growth  of  organized,  coherent  structures.3^’33  This  growth  rate  has  been  character¬ 
ized  in  semi-empirical  laws  for  incompressible  shear  layers36,37  with  corrections  for 
compressibility.37-39  If  the  aberrating  structure  had  a  characteristic  size,  A,  the  wave¬ 
front’s  frequency  of  variation,  f^pz,  would  be  expected  to  vary  with  the  structure’s 
convection  velocity  Uc  as 

/oPT=X-  (') 

The  average  coherent  structure  (a  vortex  “roller”)  diameter,  5viz,  measured  from 
schlieren  photographs  of  compressible  shear  layers,  can  be  estimated  from  the  follow¬ 
ing  empirical  relation  reported  in  the  literature:39 

f  3yiz  ^  q  (1  V)(l  ~b  VS)  ^2) 

\  x  )  5  1  +rus/J 

where  ru  is  the  velocity  ratio  U2/U1,  s  =  (h/ Pi  and  Cs  is  a  constant.  For  the  AEDC 
conditions,  Eq.  2  predicts  5viz  ~  12.6  cm  at  the  center  of  the  station  2  aperture;  that  is, 
a  structure  more  than  twice  the  size  of  the  5-cm  test  apertures  used  for  both  the  SABT 
and  the  holographic  interferometry  study  for  which  the  AEDC  facility  was  originally 
designed.33  The  effect  of  the  aperture  is  to  act  as  a  spatial  filter  on  the  measured 
wavefront.  The  distortion  caused  by  these  largest-scale  flow  structures  is  reduced  to 
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a  time- varying,  tilt  aberration  with  a  period  indicative  of  the  large-scale  structure.23,40 
The  large-scale  structures’  size  and  convection  velocity  ( Uc  ~  150  m/s)39  would  lead 
one  to  expect  tilt-distortion  variation  frequencies  of  foPZ  ~  1-2  kHz,  well  below  the 
corner  frequency  of  the  originally-used  vibration  filter.  In  the  original  study  funded 
under  our  first  grant  (and  documented  by  Ref.  18),  no  attempt  was  made  to  try 
to  understand  the  relevant  physical  cause  of  the  aberrations;  Ref.  18  sought  simply 
to  present  the  aperture-restricted  data.  The  present  grant  research,  however,  was 
specifically  interested  in  understanding  the  underlying  physics  of  compressibility  and 
anchoring  numerical  simulations  to  these  AEDC  data.  Thus  the  present  research 
required  a  more  careful  selection  of  the  high-pass  corner  frequency  in  order  to  remove 
the  vibration  corruption,  on  the  one  hand,  but  to  preserve  the  largest-scale-structure 
influences  on  the  other.  A  careful  reexamination  of  the  AEDC  data  showed  that  the 
filter’s  corner  frequency  could  be  safely  reduced  to  750  Hz  (as  described  in  Refs.  30- 
32). 

With  this  lower  vibration-removal  filter  setting,  clear  evidence  of  large-scale 
structures  was  present  in  AEDC  station  2  OPD  time-series  data  as  shown  in  Fig.  6. 
These  large-scale  flow  structures  appear  to  be  comparable  in  size  to  the  rollers  (vor¬ 
tices)  created  during  shear-layer  rollup.  For  optical  systems  with  greater  than  5-cm 
apertures,  the  structures  would  also  produce  much  larger  overall  OPDs  than  origi¬ 
nally  reported.  Since  typical  airborne  telescopes  and  laser  systems  use  apertures  much 
larger  than  the  5-cm  aperture  tested  at  AEDC,  significantly  larger  optical  degradation 
than  reported  in  Refs.  17  and  18  could  result. 

Distortions  created  by  smaller-scale  flow  structures  were  present  in  the  AEDC 
station  1  OPD  data  as  shown  in  Fig.  7.  Since  this  aperture  began  at  the  splitter-plate 
trailing  edge,  these  smaller-scale  structures  were  too  large  to  be  caused  by  shear- 
layer  rollup.  A  possible  mechanism  was  suggested  to  explain  smaller-scale  distortions 
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Figure  6.  AEDC  station  2  experimental  wavefront  reconstruction:  beam  separation 
5  =  2.54  cm. 

observed  in  both  the  station  1  and  station  2  wavefronts.  If  this  mechanism  proves 
to  be  correct,  the  overall  wavefront  might  be  modeled  as  a  superposition  of  large- 
and  small-scale  distortions  (see  below).  The  small-scale  distortions  might  then  be 
removable  by,  for  example,  boundary  layer  suction.  This  theory  is  described  in  more 
detail  in  Refs.  30-32. 

3.3.  Analytical/Numerical  Model 

To  investigate  the  physical  mechanisms  which  produce  optical  distortions  in  compress¬ 
ible  shear  layers,  optical  wavefronts  were  computed  for  a  numerical  flowfield  model. 
The  numerical  simulation’s  velocity  boundary  conditions  matched  those  of  the  AEDC 
testing  described  above.  The  simulation  had  two  distinct  parts:  (1)  a  velocity  field 
modeled  using  a  discrete  vortex  method  (DVM)  and  (2)  an  index-of-refraction  model 
used  to  determine  the  corresponding  density  field  (and  index-of-refraction  field)  from 
the  (input)  velocity  field. 
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Figure  7.  AEDC  station  1  experimental  wavefront  reconstruction:  beam  separation 
5  =  2.54  cm. 

The  initial,  transitionally-turbulent  flowfield  in  a  shear  layer  is  dominated 
by  the  shear-layer  rollup  caused  by  the  inviscid  Kelvin-Helmholtz  instability  mecha¬ 
nism.41  This  rollup  has  been  successfully  simulated  using  various  inviscid  and  pseudo- 
inviscid  methods.42-45  The  shear  layer  for  the  present  study  was  modeled  by  modifying 
the  discrete- vortex-method  (DVM)  code  developed  under  our  first  grant.  One  of  the 
primary  advantages  of  using  a  DVM  is  that  it  is  a  reduced-order,  pseudo- inviscid 
method.  As  such,  the  DVM  is  able  to  model  a  shear  flow  with  a  fairly  small  com¬ 
putational  domain,  thereby  producing  considerable  computational  savings.  A  second 
advantage  is  the  simplicity  of  the  DVM  model  itself.  This  simplicity  allows  insight  into 
the  cause-and-effect  relations  of  the  physics  being  modeled;  such  insight  was  valuable 
during  the  development  of  aero-optical  scaling  laws  discussed  in  Section  3.4  below. 
The  justification  for  using  the  (inherently  incompressible)  DVM,  its  implementation 
details,  and  the  resulting  time- varying  velocity  fields  are  presented  in  APPENDIX  D. 

Given  this  correlated  time  series  of  instantaneous  realizations  of  the  velocity 
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field  from  the  DVM,  the  problem  of  modeling  the  weakly-compressible  aberration 
mechanism  was  taken  up.  Guided  by  suggestions  in  the  literature,  density  changes  in  a 
weakly-compressible  shear  flow  were  initially  modeled  accounting  for  only  the  effect  of 
adiabatic  heating  and  cooling,  consistent  with  Morkovin’s  Strong  Reynold’s  Analogy 
(SRA). 46,47  The  SRA  suggests  that  static  pressure  fluctuations  can  be  neglected,  and 
thus  index  (i.e.,  density)  became  a  function  of  velocity  and  total  temperature  only. 
This  mechanism  in  a  total-temperature-matched  shear  layer  (like  that  at  AEDC), 
however,  produced  a  relatively  benign  index  field  and  yielded  a  time  series  of  OPDs 
that  were  totally  incompatible  with  those  measured  at  AEDC.31,40 

Forced  to  explore  other  possible  mechanisms,  a  model  was  developed  which 
accounted  for  the  pressure  gradients  which  accompany  streamline  curvature  (i.e., 
abandoning  the  SRA  altogether).  The  pressure  and  concomitant  density  fields  were 
found  by  (iteratively)  integrating  the  unsteady  Euler  equations  as  detailed  in  Ref.  40. 
Our  weakly-compressible  model  suggested  that  total  temperature  separations  akin 
to  those  produced  in  Ranque-Hilsch  tubes  might  also  be  present.  The  final  weakly- 
compressible  model  was  improved  by  incorporating  the  Ranque-Hilsch  effect.  Unlike 
the  case  of  a  compressible  boundary  layer,46  results  of  this  final  weakly-compressible 
model  showed  that  significant  pressure  fluctuations  are  present  in  a  weakly- compressible 
shear  layers,  as  illustrated  in  Fig.  8.  While  pressure  fluctuations  are  the  dominant 
factor  producing  optical  distortions,  the  temperature-separation  effect  acts  to  reduce 
the  amplitude  of  these  distortions.  Fig.  8  also  shows  the  schlierens  and  wavefronts 
produced  by  the  computed  density  fields. 

This  final  weakly-compressible  model  produced  excellent  comparisons  with 
the  large-scale,  dominant-amplitude  aberrations  observed  in  the  AEDC  experimental 
data.  Further,  the  weakly-compressible  model’s  index-of-refraction  fields  produced 
simulated  schlierens  which  matched  the  character  of  experimental  schlierens  ubiqui- 
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Figure  8.  Compressible  shear  layer  simulation  at  two  time  steps,  one  prior  to  pairing 
(left)  and  following  multiple  pairings  (right). 


tously  reported  in  the  literature  (see  Fig.  8).  Once  shear-layer  roll-up  occurred,  in 

fact,  these  large-scale  distortions  were  relatively  insensitive  to  initial  splitter-plate 
boundary-layer  thickness.31 

As  described  in  Ref.  40,  our  simulation  reasonably  models  the  large-scale  dis¬ 
tortions  measured  by  Hugo  et  al.  in  the  compressible  shear  layer18  as  shown  in  Fig.  9. 
The  small-scale  structures  measured  in  the  AEDC  OPDs  are  compatible  with  the 
notion  that  they  are  due  to  temperature  discontinuities  introduced  into  the  shear 
layer  at  the  splitter-plate  and/ or  due  to  tunnel-access- window  boundary  layers.  In 
either  case,  small-scale  distortions  would  be  superimposed  on  the  large-scale  distor¬ 
tions  predicted  by  the  weakly-compressible  model.  Evidence  of  this  is  shown  in  the 
station  1  wavefronts  measured  immediately  downstream  of  the  splitter  plate  (Fig.  9). 
Superpositioning  the  station  1  experimental  data  on  our  station  2  simulation  produces 
wavefronts  of  similar  character  to  those  measured  at  station  2  as  shown  in  Fig.  9. 

The  discrete  vortex  velocity  model  together  with  the  weakly-compressible 
mdex-of- refraction  model  described  above  produced  a  reasonable  estimate  of  the  aero- 
optical  distortions  measured  in  the  AEDC  weakly-compressible  shear  layer. 

3.4.  Resulting  Scaling  Laws 

The  DVM/weakly-compressible-model  methods,  described  above  and  anchored  to 
the  AEDC  experimental  data,  provided  a  means  of  testing  the  validity  of  candidate 
scaling  relations.  These  DVM/weakly-compressible-model  results  were  taken  as  the 

“true”  results.  The  scaling  results,  summarized  below,  are  detailed  in  Chapter  6  of 
Ref.  31. 

Aero-optical  distortions  were  found  to  scale  with  a  simple  density  ratio  for 
cases  when  the  Mach  numbers  of  the  two  constituent  flows  were  unchanged.  As 
a  consequence  of  this,  distortion  amplitudes  decreased  with  increasing  altitude,  as 


21 


Station  2  Numerical  Simulation  Station  1  Experimental  Data 


High-Pass  Filter  750  Hz 


Simulation  +  Station  1  Experimental  Data  Station  2  Experimental  Data 


Figure  9.  Comparison  of  wavefronts  from  discrete  vortex  simulation  with  superpo- 
sitioned  station  1  (experimental)  wavefront  to  reconstructed  experimental  wavefronts 
for  AEDC  facility  station  2. 
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shown  in  Fig.  10,  causing  an  increase  in  time-averaged  Strehl  Ratio,  SR.  The  scaled 
results  matched  the  “true”  results  to  within  a  fraction  of  a  percent  for  both  mean 
and  instantaneous  OPD.  Altitude  variation  changed  the  instantaneous  amplitude 
of  the  distortion  but  not  its  spatial  scales/shapes.  The  temporal  frequencies  of  the 
distortion  were  reduced  with  the  reduction  in  the  constituent  flow  velocities  required 
to  maintain  constant  Mach  with  increasing  altitude  (through  a  rescaled  time  step). 
This  suggests  that  the  correction  rate  required  by  an  adaptive-optic  system  would 
decrease  by  ~  10%  as  altitude  increased  from  13,500  ft  to  40,000  ft.31 


Figure  10.  Effect  of  altitude  on  OPDrms  (Mach  numbers  fixed)  as  computed 
by  weakly-compressible  model  (symbols)  and  using  density  ratio  scaling  law  (lines) 
(AEDC  station  2,  6^/ 2  =  8.626  mm) 

4 

An  additional  temperature  ratio  factor  was  required  to  scale  aero-optical  data 
collected  at  one  altitude  to  another  altitude  for  cases  where  the  velocity  difference  of 
the  constituent  flows  was  held  constant.  Data  scaled  with  this  relation  were  within  4% 
of  the  results  obtained  for  the  “true”  results  both  for  instantaneous  and  mean  OPDs. 
Distortion  variation  with  altitude  followed  the  same  trends  as  for  the  constant-Mach 
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cases  as  shown  in  Fig.  11;  however,  distortion  amplitudes  were  larger  for  the  constant- 
velocity-difference  cases  due  to  the  increase  in  Mach  numbers  which  accompanies  the 
reduction  in  the  sonic  speed  with  increasing  altitude.  In  this  case,  wavefront  temporal 
frequencies  were  unaffected  by  changes  in  altitude.31 


Figure  11.  Effect  of  altitude  on  OPDrms  (velocity  difference  fixed)  as  computed  by 
weakly-compressible  model  (symbols)  and  using  combined  density/temperature  ratio 
scaling  law  (lines)  (AEDC  station  2,  Si/2  =  8.626  mm) 

The  OPDTrns  data  for  the  weakly-compressible  model  at  different  altitude  con¬ 
ditions  were  compared  to  free-shear-layer  (fence)  and  turbulent-boundary-layer  data 
obtained  in  actual  flight  tests  by  Gilbert.48  Incredibly,  the  weakly-compressible-model 
data  agree  with  Gilbert’s  data  to  within  the  experimental  measurement  uncertainty 
as  shown  in  Fig.  12.  In  addition,  a  quadratic  fit  of  the  weakly-compressible-model 
data  provided  a  superior  match  over  the  linear  relation  suggested  by  Gilbert,  although 
the  general  trends  of  both  fits  are  similar.31 

The  scaling  objective  of  this  grant  was  prompted  by  the  inability  to  scale  the 
OPD  structure  sizes/frequencies  of  the  heated  jet  measurements  to  those  measured 
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Figure  12.  Effect  of  dynamic  pressure,  on  OPDrrns.  KC-135  data  taken 
from  Ref.  48;  DVM/weakly-compressible-model  data  from  AEDC  station  2,  5l/2  = 
8.626  mm  simulation. 

at  AEDC.  The  present  grant’s  research  showed  that  AEDC  flowfield  structure  size, 
dominant  temporal  frequency,  and  maximum  distortion  amplitudes  scaled  as  would 
be  expected  from  the  fluid-mechanics  literature.  This  scaling  is  based  on  the  primary 
coherent  structures  in  the  shear  layer  and  thus  only  holds  for  an  aperture  sized  larger 
than  the  largest  flow  structures.  Since  the  AEDC  data  was  obtained  with  an  aperture 
smaller  than  this  structure  size,  the  scaling  was  only  apparent  by  comparing  the 
AEDC- anchored  DVM/weakly-compressible-model  results  with  those  of  the  heated 
jet  as  shown  in  Table  1.  The  specifics  of  this  scaling  are  detailed  in  Section  6.5  of 
Ref.  31. 


4.  TWO-DIMENSIONAL  WAVEFRONT  SENSOR 

In  the  first  grant,  we  demonstrated  a  high-bandwidth  100  kHz)  sensor  which  took 
advantage  of  the  fact  that  aberrating,  turbulent  flow  structures  “convect”  across 
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Table  1.  Comparison  of  flowfield  characteristics  for  two  example  shear  flows:  heated 
jet  (dissimilar- index  mixing)  and  weakly-compressible  shear  layer  (DVM/weakly-com- 
pressible  model) 


Facility 

Heated  Jet 

AEDC  Shear  Layer 

Ratio 

(Ref.  6) 

=  AEDC/ Jet 

Distortion  Cause 

M\  ~  0.8, 

WM 

o 

l 

£ 

Tested  Aperture  (cm) 

1.27 

20 

~  16 

Equivalent  D  (cm) 

1.27 

20 

~  16 

Equivalent  x  (cm) 

3.18 

48.3 

~  15 

Equivalent  x/D 

2.5 

2.4 

~  1 

AU  (m/s) 

7 

226 

-  30 

Ue  (m/s) 

3.5 

148 

~  40 

Structure  Size,  A  (mm) 

6 

O 

o 

T— 1 

l 

~  16 

foPL  (Hz) 

~  500 

~  1500 

~  3 

OP Dmax  (^m) 

0.1 

~  1.5 

~  15 

the  viewing  aperture.6,13,18  This  reliance  on  structure  convection,  however,  limits  the 
SABT  to  measurement  of  the  streamwise  component  of  the  optical  wavefront.  For  the 
two-dimensional  derivative  of  the  SABT,  we  used  the  streamwise  convection  to  unfold 
the  wavefront  in  the  streamwise  direction.  In  addition,  we  imposed  a  scan  velocity 
in  the  cross-stream  direction  to  unfold  the  remaining  orthogonal  component,  thereby 
producing  the  full,  two-dimensional  wavefront  as  shown  schematically  in  Fig.  13. 

The  testing  of  a  breadboard  implementation  of  the  cross-stream  scanner  com¬ 
ponent  revealed  a  measurable  beam  displacement  inherent  in  the  scanning  process 
alone.14  This  displacement  can  be  removed  as  a  tare,  provided  it  is  repeatable.  Thus, 
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(a)  (b) 


Figure  13.  Example  (simulation)  of  (a)  a  two-dimensional  flowfield  and  (b)  the 
resulting  two-dimensional  optical  wavefront  over  the  aperture  defined  by  the  “up¬ 
stream”  and  “downstream”  beams.  (From  Ref.  15.) 

a  key  to  successful  implementation  of  this  instrument  concept  was  the  incorpora¬ 
tion  of  a  galvanometer  scanner  with  position  feedback  control  for  high  scan-to-scan 
repeatability.15 

The  SABT-derivative  instrument  was  tested  in  the  Notre  Dame  water  tunnel. 
A  mixing  layer  was  produced  by  a  vertically-mounted  splitter  plate  at  the  upstream 
end  of  the  test  section  as  shown  in  Fig.  14  and  detailed  in  Ref.  15.  A  velocity 
differential  between  the  two  sides  of  the  plate  was  created  by  installing  a  “straw 
box”  flow  restrictor  on  one  side  of  the  plate  and  by  modifying  the  tunnel  inlet.  A 
1000-W  cartridge  heater,  installed  11.4  cm  ahead  of  the  trailing  edge,  produced  an 
optically-active  flow.  During  this  experiment,  the  flow  on  one  side  of  the  plate  was 
Ui  =  18.5  cm/sec  while  the  flow  on  the  low-speed  side  was  U2  =  4  cm/sec.  The 
flowfield  is  visualized  with  hydrogen  bubbles  in  Fig.  15.  Visible  in  Fig.  15  are  vortical 
coherent  structures  in  the  mixing  layer  itself  which  are  roughly  equivalent  to  those 
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shown  the  two-dimensional  flow  example,  Fig.  13.  The  resulting  wavefronts  were 
a  time-resolved,  time  series;  individual  wavefronts  were  obtained  at  600  Hz  (two 
measurements  per  scanner  cycle).15  A  series  of  these  wavefronts  (every  20th  frame) 
is  shown  in  Fig.  16.  The  bandwidth  of  this  new  wavefront  sensor  can  be  increased  to 
8  kHz  by  incorporating  a  faster  scanner  (already  in  hand).  We  plan  to  use  the  8  kHz 
instrument  for  two-dimensional  wavefront  measurements  in  laboratory  air  flows. 


Y 


Water  Tunnel 
Test  Section 


(a) 


(b) 


Figure  14.  (a)  Photograph  of  water  tunnel  heated  shear  layer  facility  and  (b) 

schematic  of  the  SABT-derivative  instrument  optical  bench. 


5.  SHEAR-LAYER  CONTROL 

The  use  of  control  for  heated-jet  experiments  greatly  enhanced  certain  areas  of  in¬ 
vestigation  in  our  previous  work  (see,  for  example,  Refs.  49  and  50).  The  coupling  of 
control  with  the  SABT  sensor  suggests  the  possibility  of  not  only  real-time  sensing 
of  wavefront  aberrations  but  also  for  adaptive-optic  correction.  An  initial  look  at 
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Figure  15.  Hydrogen  bubble  visualization  of  heated  mixing  layer  (looking  up 
through  the  bottom  of  tunnel  test  section). 
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Figure  16*  Experimental,  two-dimensional,  wavefront  reconstructions:  time  steps 
=  924.  (a)  through  1024  (f).  Between  frames,  At  =  0.0333  sec;  overall  period  (a)-(f) 


0.1667  sec. 
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the  possibility  of  shear-layer  control  was  performed  during  this  grant  using  the  DVM 
code. 


For  a  free  shear  layer,  the  wavelength  of  the  most  unstable  mode  is  driven 
by  the  layer’s  initial  thickness  according  to  linear  stability  theory51’52  and  the  exper¬ 
imental  literature.53,54  In  fact,  we  found  that  our  DVM  model  followed  the  theory 
extremely  well.  The  shear  layer’s  natural  frequency,  can  be  predicted  from  the 
Strouhal  number 

St0  —  yy  —  0.032  (3) 

where  90  is  the  momentum  thickness  of  the  high-speed  stream  at  the  splitter  plate 
trailing  edge  (and  assuming  Uc  =  (Ui  +  f/2)/2).53’54  Shear  layer  rollup  has  been  exper¬ 
imentally  observed  to  occur  at  approximately  x  =  3An  where  the  natural  wavelength, 
An,  is  given  by53’54 

A„  =  y  •  (4) 

In 

Table  2  compares  the  shear  layer  natural  frequencies  and  concomitant  first  rollups 
predicted  using  these  relations  for  different  initial  discrete  vortex  core  radii.  (It  should 
be  noted  that  the  length  scale  0o  used  in  Table  2  was  the  momentum  thickness  of  the 
high-speed  stream  only,  as  stipulated  in  Ref.  53.)  Table  2  shows  that  the  DVM’s  initial 
rollup  occurs  at  a  position  that  varies  with  initial  layer  thickness  and  is  consistent 
with  theory  and  experiment.  Moreover,  if  the  shear  layer  is  artificially  perturbed  at 
frequency  /„,  the  shear  layer  tends  to  rollup  into  “rollers”  of  length  ~  An  as  shown 
in  Fig.  17. 

The  streamwise  energy  content  of  the  forced  shear  layer  is  typically  quantified 


as' 


53,54 


K(/)F 

20o 


dy 


(5) 


where  u'(f )  is  the  narrow-band  velocity  fluctuation  level.  When  the  energy  content 
within  each  frequency  band  was  computed  from  forced  DVM  velocity  profiles,  the 
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Si 

(mm) 

#0 

(mm) 

1 

A„ 

(cm) 

Predicted 

Rollup  (cm) 

DVM  Typical 

Rollup  (cm) 

69.00 

5.250 

901.29 

16.41 

49.22 

-  45-60 

34.50 

2.975 

1590.80 

9.30 

27.89 

~  26-30 

17.25 

1.429 

3310.79 

4.46 

13.40 

~  16-20 

Table  2.  DVM  simulation  of  (unforced)  AEDC-type  shear  layer:  DVM  rollup  posi¬ 
tion  compared  to  predictions  using  Ref.54  methods  (0O  for  high-speed  stream  only) 

energy  at  each  frequency  initially  increased  with  increasing  x ,  as  shown  in  Fig.  18 
for  frequencies  /„  and  /„/ 2.  The  energy  content  in  the  forcing  frequency  mode,  /„, 
increased  until  saturation  occurred  at  the  x  corresponding  to  the  shear-layer  initial 
rollup.  Beyond  the  initial  rollup,  the  energy  content  at  the  subharmonic  frequency 
/„/ 2  dominated  until  it  also  saturated  at  the  position  corresponding  to  the  first  pairing 
of  the  vortex  “rollers.”  This  energy  amplification  behavior  matches  that  noted  in  the 
experimental  literature. 


Figure  IT.  Locus  of  discrete  vortex  positions  for  artificially  forced  DVM  shear  layer, 
time  step  =  716  (AEDC  conditions,  Si  =  34.5  mm,  forcing  frequency,  fn  =  1591  Hz) 

The  good  agreement  between  the  DVM  simulation  and  the  behavior  of  real 
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Figure  18.  Streamwise  energy  content  variation  with  x  for  artificially  forced  DVM 
shear  layer  (AEDC  conditions,  Si  =  34.5  mm,  forcing  frequency,  /„  =  1591  Hz) 

shear  layers  reported  in  the  experimental  literature  suggest  that  the  DVM  is  a  useful 
tool  for  evaluating  potential  forcing  schemes. 

The  problem  of  flow  control  at  compressible  velocities  centers  on  having  suf¬ 
ficient  frequency  response  and  actuator  energy  density.  Electro-menchanical  devices 
which  can  produce  high  enough  amplitudes  (electro-magnetic  or  bi-metallic  pistons 
and  flaps,  for  instance)  typically  lack  sufficient  bandwidth  whereas  those  with  high 
frequency  response  (piezoelectric,  electrostatic)  typically  have  low  displacement  am¬ 
plitudes.  Corke  and  Cavalieri55,56  have  used  phased-plasma-array  actuators  for  excit¬ 
ing  instabilities  in  a  Mach  3.5  boundary  layer.  These  show  promise  for  shear-layer 
control  as  well. 

This  problem  of  shear-layer  control  will  be  further  explored  both  numerically 
and  experimentally  in  our  follow-on  grant,  AFOSR  F49620-00- 1-0025. 


32 


6.  INFERRING  FLUID-MECHANIC  INFORMATION 

FROM  WAVEFRONTS 

Finally,  some  preliminary  investigations  into  inferring  fluid-mechanic  information 
from  the  distorted  wavefronts  were  undertaken  in  our  first  grant  (see,  for  example, 
Refs.  21  and  57);  however,  in  this  objective,  we  have  begun  to  investigate  this  rela¬ 
tionship  more  formally  by  employing  wavelet  transforms  and/or  proper  orthogonal 
decomposition  methods  to  decompose  the  SABT  signals.  In  so  doing,  we  plan  to  de¬ 
velop  algorithms  for  extracting  the  velocities  for  each  of  the  aberration-causing  flow 
structures.  This  work  is  the  primary  focus  of  our  newest  researcher  and  will  continue 
in  our  follow-on  grant,  AFOSR  F49620-00- 1-0025. 

7.  DISSEMINATION  OF  RESULTS 

7.1.  Publications 

This  research  has  produced  a  Ph.D.  dissertation  (Ref.  31),  six  meeting  papers  (Refs.  21, 
14,  15,  58,  40,  and  30),  and  two  journal  articles  (Refs.  16,  59).  Refs.  60  and  32,  versions 
of  Refs.  15  and  30  (meeting  papers),  are  currently  in  review  by  the  journals  Optical 
Engineering  and  AIAA  Journal ,  respectively.  Ref.  40  is  currently  being  revised  for 
submission  to  Journal  of  Fluid  Mechanics.  In  addition,  a  tentative  arrangement  has 
been  reached  with  the  editor  of  Progress  in  Aerospace  Sciences  to  publish  an  updated 
version  of  the  “AIAA  Plasmadynamics  and  Lasers  Best  Paper  Award  (1997)”-winning 
review  paper  (Ref.  21)  funded  under  this  grant. 

7.2.  Technology  Transfers 

In  addition  to  presentations  made  at  AFOSR-sponsored  meetings  and  national  tech¬ 
nical  meetings,  we  have  endeavored  to  transition  our  findings  to  those  who  are 
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actively  working  in  the  aero-optics  field.  Recently,  our  research  was  presented  by 
Dr.  D.  Kyrazis  at  the  Second  Annual  Directed  Energy  Symposium  (1-5  November 
1999). 19  In  addition,  we  have  explored  with  Boeing  the  use  of  the  SABT  for  ABL 
aero-optics  ground  tests  and  have  offered  to  assist  in  developing  an  array  of  sensors 
for  this  purpose.  We  have  also  supported  the  hypersonic  aero-optic  testing  performed 
in  Calspan  UB  Research  Center’s  Large-Energy  National  Shock  facility. 

8.  PERSONNEL  SUPPORTED  UNDER  THIS  EFFORT 

This  grant  provided  approximately  one  month  of  support  each  year  to  the  Principal 
Investigator,  Dr.  Eric  J.  Jumper.  In  addition,  a  full-time  Ph.D.  student,  Edward  J. 
Fitzgerald,  was  supported  through  his  Ph.D.  degree.  A  second  Ph.D.  student,  James 
Cicchiello,  was  also  supported  during  the  last  year  of  the  grant. 
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APPENDIX  A.  AERO-OPTIC  STATISTICAL  APPROACH 

An  early,  if  not  the  first  attempt  to  quantify  the  aero-optical  aberration  imposed  by 
propagation  of  a  collimated  coherent  beam  of  light  through  an  optically- active  (index- 
of-refraction  variant)  turbulent  flowfield  was  that  by  Liepmann.8  The  motivation  for 
this  first  study  was  to  place  some  limits  on  the  sensitivity  (i.e.,  sharpness)  of  Schlieren 
systems  for  use  in  high-speed  flow  analysis.  Liepmann  analyzed  the  propagation  of  a 
small-aperture  (small-diameter)  beam  and  derived  its  mean-square  deflection  angle  off 
its  mean  angle,  (e2),  using  geometric  optics.  After  propagating  through  a  boundary 
layer  of  thickness  d,  in  the  propagation  direction  y,  where  changes  in  the  index  of 
refraction  in  the  normal  plane,  x-z,  are  taken  to  be  statistically  similar,  Liepmann’s 
analysis  gave 

rfc  (6) 

where  the  index  of  refraction  is  given  by  n  =  n0(t/)(l  +  v ),  and  Rv(\  y-(  |)  is  the  cor¬ 
relation  function  for  the  index  variation.  When  the  index-of  refraction  fluctuation,  v, 
is  replaced  by  the  density  fluctuation,  p',  using  the  Gladstone-Dale  constant,  Kqd ,22 

p'Kgd  ( *7\ 

v= -  (7) 

n0 


Eq.  6  becomes 


Kgd 

n0(5) 


f  f( 

Jo  Jo 


)  RP' (I  y  -C I)  dv  dC 


It  is  of  interest  to  note  that  this  early  work  introduced  concepts  that  remain  important 
to  aero-optics;  first,  Eqs.  6  and  8  introduce  the  notion  of  linking  statistical  measures 
of  the  turbulent  flowfield  to  the  time-averaged  optical  degradation.  Further,  they 
introduce  the  correlation  function  or  “scale  of  the  turbulence”  and  its  position  depen¬ 
dence  along  the  propagation  direction  as  key  to  evaluating  the  statistical  aero-optical 


effect. 
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A  significant  piece  of  work  was  done  in  1956  by  Stine  and  Winovich;61  they 
performed  photometric  measurements  of  the  time-averaged  radiation  field  at  the  focal 
plane  of  a  receiving  telescope  that  imaged  an  initially  planar  light  wave  apertured  by 
pinholes  ranging  in  diameter  from  0.75  x  10"3  to  33.0  x  10-3  inches.  This  experimental 
data  was  used  to  assess  Liepmann’s  analysis  as  well  as  that  performed  at  Northrop 
Aircraft  by  Baskins  and  Hamilton  on  the  transmission  characteristics  of  an  optical 
signal  through  supersonic,  turbulent  boundary  layers.62,63  Stine  and  Winovich’s  work 
is  important  on  several  accounts.  It  brought  together  all  that  had  been  done  till 
then  on  optical  propagation  through  index- variant  turbulent  flow.  In  addition,  it 
brought  to  bear  on  the  problem  the  relatively  mature  analysis  of  electro-magnetic 
scattering  in  the  troposphere  and  through  atmospheric  turbulence64,65  which  would 
continue  to  be  of  importance  to  aero-optics.  This  work  also  raised  the  prospect 
of  using  an  optical  degradation  measurement  as  a  method  of  inferring  turbulence 
scales;  this  hope  of  being  able  to  infer  information  about  the  turbulent  flowfield 
from  optical  measurements  of  aero-optic  effects  remains  as  one  of  the  potential  fluid- 
dynamic  payoffs  that  may  come  from  aero-optic  research. 

In  my  opinion  the  milestone  paper  in  “modern”  aero-optics  is  that  by  Sutton;2 
this  paper  set  several  trends  that  dominated  the  study  of  statistical  aero-optics  in  the 
“modern”  era.  By  “modern”  I  mean  here  the  effect  of  aero-optics  on  the  beam  quality 
of  high-power  laser  irradiation.  The  aero-optic  problems  in  the  laser  applications  are 
ubiquitous,  and  the  combined  fluid-mechanist  and  laser-physicist  communities  began 
to  get  serious  about  dealing  with  the  problems.  The  Sutton  paper  was  the  first 
paper  to  consider  the  problem  of  the  effect  of  near-field  phase  aberrations  due  to 
the  propagation  of  a  large-aperture,  single  wavelength,  collimated  (planar-wavefront) 
beam  through  a  turbulent,  optically-active,  spatially-confined  flowfield.  Imbedded  in 
the  analysis  was  the  first  time  the  linking  equation  relating  optical  phase  variance 
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1  • 


over  the  exit  pupil  to  specific  fluid-mechanic  properties  was  advanced  in  the  form  of 
a  position  dependent  (along  the  propagation  path)  extinction  coefficient,  a, 

a  =  2 k2  ( v 2)  A  (9) 


in  the  Lambert’s  law  sense  for  a  reduction  in  intensity,  7,  as 


fy  -  ~aI  <10> 

where  k  =  27r/A,  A  being  the  wavelength  of  the  radiation.  When  integrated  over  the 
thickness  of  the  flowfield,  L,  Eq.  9  yielded  the  mean-square  phase  error,  as 

al  =  2k2J0  al'lvdy  (ii) 

where  Sutton  suggested  that  the  proper  correlation  length  scale,  ly)  should  be  A, 
the  position-dependent  integral  scale  size  of  the  turbulence.66  If  the  Gladstone-Dale 
constant  is  introduced  as  above,  Eq.  11  becomes 

oj  =  2 K2gd  k2  Jo  °2pf  My)  dy •  (12) 

The  paper  also  introduced  the  now  common  use  of  the  mean-square  phase  error  to 
directly  estimate  of  the  far-field  intensity  reduction  or  Strehl  ratio1 


SR  = 


/(Q,o) 

/o(0,0) 


(13) 


via  the  large-aperture  approximation  to  obtain  the  time-averaged  Strehl  ratio  as 


SR  =  exp(-aj). 


(14) 


From  the  time  of  that  paper  on,  the  focus  and  content  of  aero-optic  research 
was  shaped.  A  review  of  the  major  publication  in  aero-optics,  Aero- Optical  Phe¬ 
nomena ,67  demonstrates  that  work  up  until  1982  focused  on  the  measurement  of  the 
time-averaged,  spatial,  near-field  optical  phase  variance,  Measurements  of  were 
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obtained  either  by  direct  optically-based  methods  (Ref.  68,  for  example),  or  indirectly 
using  fluid-mechanic  measurements  via  the  linking  equation  of  Eq.  12  (Ref.  69,  for  ex¬ 
ample).  Optical  methods  that  have  been  applied  to  the  measurement  of  the  near-field 
time-averaged  phase  variance  include  direct  interferometry,  pulsed  interferometry68 
and  shearing  interferometry.70  These  interferometric  methods  provide  a  time-averaged 
assessment  of  the  optical  phase  variance  over  the  aperture  by  quantifying  the  phase 
variance  of  individual  interferograms  and  averaging  the  variance  over  the  number  of 
interferograms  used  to  quantify  the  aperture  aberrations.  Individual  interferograms 
are  uncorrelated,  so  that,  given  a  sufficient  number  of  interferograms,  the  estimates  of 
the  time-averaged  phase  variance  can  be  expected  to  yield  accurate  results;  however, 
these  methods  provide  no  information  concerning  temporal  frequencies. 

Indirect  methods  of  assessing  the  near-field  optical  phase  variance  via  Eq.  12 
depend  on  a  quantification  of  the  density  fluctuations  which  are  generally  experi¬ 
mentally  measured  using  hot-wire  methods.2,69  Since  a  complete  quantification  of  a 
two-point  correlation  is  time  consuming,  Rose  et  alm  developed  further  assumptions 
regarding  the  form  of  the  covariance  function  of  the  density  fluctuations  to  reduce 
the  required  sampling  time  of  the  measurements.  Again  these  methods  provide  no 
temporal  information  about  the  near-field  phase  variance.  Many  efforts  were  directed 
at  comparing  the  inferred  phase  variance  from  fluid-mechanic  linking  equations  and 
the  optically-measured  (inferred)  time-averaged  phase  variance.  As  recently  as  1992, 
a  paper  was  presented  reviewing  the  analytic  linking  methods  used  to  infer  optical 
wavefront  variance.9 

To  end  this  brief  and  incomplete  historical  review,  I  refer  to  a  final  paper, 
again  by  Sutton,  published  in  1985.66  I  do  not  presume  to  judge  the  intent  of  that 
paper;  however,  to  this  reader,  the  paper  suggests  that  the  field  of  aero-optics  was 
essentially  mature,  and  that  the  opportunities  for  further  research  were  to  polish  up 
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the  techniques  already  in  existence  and  maybe  apply  them  to  a  few  more  flows.  I 
should  note  that,  in  fact,  work  continues  to  this  day  on  dotting  the  i’s  and  crossing 
the  t’s  of  the  traditional  methods  of  assessing  near-field  phase  variance.  In  the  end, 
however,  the  now-traditional  methods  in  aero-optics  are  limited  in  their  ability  to 
take  an  active  role  in  system  design  and  performance.  As  the  field  of  aero-optics  has 
been  described,  the  best  that  can  be  hoped  for  is  an  ever  more  accurate  assessment  of 
the  optical  degradation  that  may  be  expected  (in  a  time-averaged  sense)  for  a  given 
set  of  flow  conditions.  These  methods  do  not  suggest  either  the  temporal  or  spatial 
frequencies  associated  with  the  optical  aberrations  and,  more  importantly,  they  do  not 
even  provide  a  close  approximation  of  the  instantaneous  near-field  optical  wavefront 
distortions.  In  terms  of  providing  design  specifications  for  adaptive-optical  systems, 
the  traditional  aero-optics  approaches  are  of  little  to  no  value;  in  other  words,  as 
represented  up  to  1992,  aero-optics  was  at  a  level-one  understanding  and  content  to 
be  there. 
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APPENDIX  B.  FIRST  DIRECT  OPTICAL 
MEASUREMENT  OF  OPDrms 

The  first  instrument  to  address  dynamics  of  aero-optics  was  developed  by  Malley  et  al. 
in  the  late  1980’s,10  although  their  intention  was  only  to  obtain  a  statistical  measure  of 
OPDrms.  Their  concept  first  appeared  in  print  in  1989.11  In  these  papers,  Malley,  et 
al.  describe  the  invention  of  a  relatively  inexpensive  and  quick  optical  method  for  as¬ 
sessing  the  statistical,  optical  phase  variance.  Key  to  their  method  was  Malley  et  al.'1  s 
assertion  that  optical  aberrations  imposed  by  a  convecting,  optically-active  turbulent 
flow  must  themselves  “convect”  through  the  viewing  aperture.  In  their  writings  this 
assertion  is  stated  almost  cavalierly,  as  if  it  were  common  knowledge,  and  so,  upon 
thinking  about  it,  it  would  seem;  however,  this  assertion  must  be  cited  as  one  of  the 
pinnacle  discoveries  of  the  new  dynamic  era  in  aero-optics.  With  this  recognition, 
a  spatially-stationary  time  series  of  wavefront,  spatial-derivative  measurements  can 
be  used  to  determine  an  accurate  statistical  measure  of  the  wavefront  distortion  over 
a  relatively-large  region  of  the  viewing  aperture  near  the  point  where  the  derivative 
measurement  is  made.  I  will  return  specifically  to  Malley  et  aVs  method;  however, 
first  I  want  to  explore  the  ramifications  of  their  discovery. 

B.l.  Hartmann  Sensor  Theory 

A  common  method  of  measuring  the  figure  of  an  optical  wavefront  is  by  employing  a 
Hartmann  wavefront  sensor.12  According  to  Huygen’s  principle71  a  ray  of  light  travels 
normal  to  its  associated  optical  wavefront;  this  principle  is,  in  fact,  the  basis  for 
geometric  optics.  As  shown  in  Fig.  19  (greatly  exaggerated),  if  a  planar  wavefront 
is  made  to  propagate  through  an  aberrating  medium  (due  to  density  fluctuations, 
for  example),  the  wavefront  emerges  distorted  or  aberrated.  Note  in  Fig.  19  that 
there  are  a  number  of  vectors  shown  normal  to  the  emerging,  aberrated  wavefront, 
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Figure  19.  Planar  wavefront  and  probe  beam  distorted  by  a  variant  fluid  density 
field  (left)  and  resulting  Hartmann  Sensor  dot  displacements  (right). 


and  also  the  path  of  a  single  “ray”  as  it  traverses  the  medium  and  emerges  normal 
the  wavefront.  In  a  classical  Hartmann  Plate,  an  opaque  plate  with  an  array  of 
small-diameter  (small-aperture)  holes  in  it  is  placed  at  the  point  where  the  wavefront 
emerges  from  the  aberrating  medium,  forming  as  many  small-aperture  beams  as  there 
are  holes  in  the  plate.  Each  hole  in  the  Hartmann  Plate  defines  the  centroid  of  a 
subaperture  of  the  full  viewing  aperture.  Because  of  Huygen’s  principle,  these  beams 
will  emerge  from  the  plate  at  an  angle  normal  to  the  plate-impinging  wavefront  at 
the  location  of  the  hole.  As  shown  in  Fig.  19,  If  a  photographic  plate  is  placed 
at  some  distance  d  from  the  plate  in  the  y  direction,  and  exposed  for  a  very  short 
time,  the  photographic  plate  will  record  the  displacement  of  each  beam  from  its  null 
position.  The  null  position  can  also  be  recorded  onto  the  photographic  plate  by 
exposing  the  plate  to  the  small-aperture  beams  emerging  from  the  Hartmann  Plate 
with  no  aberrating  medium  present.  An  example  of  the  unaberrated  and  aberrated 
photographic  plate  might  look  like  those  in  Fig.  19.  If  the  distance  in  the  y-direction 
of  the  photographic  plate  from  the  Hartmann  plate  is  d,  and  the  displacement  of  the 
dot  from  one  small- aperture  beam  from  its  null  position  in  the  x  direction  is  ex,  then 
6X.  the  off-planar  angle  of  the  wavefront  in  the  x-y  plane,  at  the  location  where  the 
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wavefront  impacts  the  Hartmann  plate  can  be  determined,  referring  to  Fig.  19,  by 

6X  =  arctan(^).  (15) 

If  9X  is  small,  as  it  will  normally  be,  then 

0*  =  (16) 

a 

9X  is  the  derivative  of  the  wavefront  in  the  x-direction  in  the  x-y  plane  normal  to 
the  incoming  wavefront  direction  at  the  Hartmann  Plate.  As  described  in  Ref.  12, 
the  relative  wavefront  (relative  Optical  Path  Length,  OPL)  in  the  x  direction  can  be 
constructed  approximately  by 

0PL(Xj)  ~  1 E  (£xn~12+  c*n')  ~  Xn~i)’  (1?) 

where  xx  is  the  location  of  the  first  hole  in  the  Hartmann  Plate  at  x  =  0.0,  x2  in  the 
next  location,  etc.,  and  Xj  is  the  location  of  the  jth  of  N  holes  in  the  x-direction  in  the 
Hartmann  Plate.  The  reason  that  this  equation  is  approximate  is  that  it  represents 
a  discretized  approximation  to  the  integral 

(18) 

The  validity  of  the  approximation  depends  on  the  distance  between  holes,  (xn  —  xn_i), 
being  small  compared  to  the  change  in  the  slopes  of  the  wavefront.  A  similar  approach 
using  the  displacement  of  the  dots  on  the  photographic  plate  from  their  null  positions 
in  the  z-direction,  ez,  gives  OPL(z).  Complete  constructions  of  OPL(x)  and  OPL(z) 
over  the  viewing  aperture  yields  a  two-dimensional  wavefront  at  the  instant  in  time 
at  which  the  photographic  plate  is  exposed,  OPL(t,x,  z). 

If,  instead  of  a  photographic  plate,  the  plate  is  replaced  with  a  charge-coupled 
device  (CCD)  array,  and  the  Hartmann  Plate  is  replaced  by  an  array  of  lenslets  (each 
lenslet  defining  a  subaperture),  the  wavefront  sensor  is  termed  a  Shack-Hartmann 
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sensor.  This  sensor  is  capable  of  making  wavefront  measurements  at  the  framing 
rate  of  the  CCD  array.  The  standard  video  framing  rate  is  30  Hz,  but  research- 
class,  state-of-the-art  CCD  arrays  can  frame,  i.e.,  interrogate  every  pixel  in  the  CCD 
array,  in  approximately  1/1000  sec,  or  a  rate  of  1  kHz,  to  produce  a  time  series 
of  OPL(x,  z,t).  Such  wavefront  sensors  are  now  in  common  use  in  adaptive  optic 
systems  for  the  atmospheric  propagation  problem. 


Malley  et  a/.’s  discovery  recognized  that  if  the  aberrations  were  convecting 
(unchanging),  then  the  required  closely-spaced  individual  measurements  of  the  wave- 
front  derivative,  in  the  form  of  off-axis,  small-aperture-beam  displacements  could  be 
replaced  by  a  time  series  of  off-axis  displacements  at  a  single  location.  This  time 
series  could  then  yield  an  optical  pathlength  as  a  function  of  time  because 


where 


is  just  —  0x(t),  the  continuous  time  series  of  off-axis  displacement  angles 


(given  by  Eq.  16)  in  the  convection  direction,  and  Uc  is  the  ^  of  the  “convecting” 
wavefront  aberration.  The  discretized  version  of  Eq.  19  is 

N 

QPL(tN) «  £  -  [ex(u)uc  At] .  (20; 


Given  the  relation  of  Eq.  19,  Hugo22  reexamined  Liepmann’s  linking  equation 
of  Eq.  8  (of  Appendix  A).  Using  a  Taylor  hypothesis  so  that  the  derivative  in  Eq.  8, 


for  example,  can  be  replaced  by 


8J_  _  _J ._dj_ 
dz  Uc  dt 


and  incorporating  Malley  et  al.’s  discovery  relating  local  tilt  and  Uc  and  integrating 
with  time  to  arrive  at  OPL(t ),  an  equation  can  be  developed  following  Liepmann’s 


derivation  of  the  form 


kgd  k2  f  [  Co v(y  -  C)  dy  dC- 

Jo  Jo 
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When  the  definition  of  the  covariance  function  is  inserted  into  Eq.  22  and  integrated , 
it  becomes 

=  2Kgd  k2  JQ  A(y)  dV  (23) 

which  is  seen  to  be  identical  to  Eq.  12  of  Appendix  A.  Thus,  Malley  et  aV s  discovery 
had  implications  that  “unified”  the  traditional  approaches  to  aero-optics  as  well  as 
marking  a  turning  point  in  the  study  of  aero-optic  phenomena. 

B.2.  The  Malley  et  al.  Device. 

In  Ref.  10  a  device  developed  by  Malley  et  al.  was  described  that  was  able  to  make 
direct  optical  measurements  of  optical  phase  variance  in  near  real  time.  As  mentioned 
above,  the  device  was  based  on  their  assertion/discovery  that  the  optical  aberrations 
caused  by  an  optically-active  turbulent  fluid  medium  convect  with  the  convecting  flow 
structures.  The  device  was  cleverly  simple  and  made  use  of  an  off-the-self  electronic 
autocollimator.  By  monitoring  the  voltage  fluctuations  from  the  collimator  they 
were  able  to  record  a  time  series  of  the  off-axis  displacement  (jitter)  of  the  centroid 
of  a  small  aperture  laser  beam.  The  beam  originated  at  the  autocollimator  and 
returned  from  a  reflection  off  a  flat  mirror  on  the  opposite  side  of  the  aero-optic 
disturbance.  As  discussed  above,  this  beam  jitter  was  directly  related  to  beam  off- 
axis  angle  variations  ( 9X  and/or  9Z).  This  small  angle  is  the  spatial  derivative  of  the 
wavefront  at  the  location  of  the  small-aperture  probe  beam  at  the  time  associated 
with  that  measurement  in  the  time  series.  As  pointed  out  by  Malley  et  al.,  once 
the  time  series  is  available,  Fourier  techniques  provide  a  direct  measure  of  the  OPD 
variance  as  a  function  of  frequency. 

Invoking  a  Taylor  frozen-flow  hypothesis  allows  for  the  spatial  frequency,  k, 
to  be  obtained  from  the  temporal  frequency  via  k  =  u>/Uc,  where  u  is  the  temporal 
frequency.  Thus,  the  instrument  provided  a  measure  of  OPD  variance  both  as  an 
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integrated  total  and  as  a  function  of  spatial  frequency. 

Although  statistical  (or  perhaps  better  stated  as,  on  a  statistical  basis),  the 
data  obtainable  from  the  Malley  et  al.  device,  contained  far  more  aero-optical  infor¬ 
mation  than  had  ever  been  previously  available.  It  also  compared  well  (was  validated 
by  comparisons)  with  the  direct  optical  measurement  of  the  OP Drms  inferred  from 
averaging  many  instantaneous  realizations  of  the  flowfield  via  interferometry.  In  one 
step,  aero-optics  had  advanced  to  a  level-two  understanding  and  completely  bypassed 
the  requirement  for  a  linking  equation  like  Eq.  12.  Later  work  by  Hugo22  confirmed 
the  high  accuracy  of  the  Malley  et  al.  device  in  measuring  OPDrms. 
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APPENDIX  C.  THE  SABT  SENSOR  DEVELOPED 
UNDER  THE  FIRST  FUNDING  ROUND 

A  new,  very  high-bandwidth  wavefront  sensor  was  developed  under  our  first  round 
of  AFOSR  funding.  Refs.  6  and  13,  describe  a  technique  of  making  use  of  Mal- 
ley  et  a/.’s  convecting  aberration  discovery,  and  multiple  probe  beams  to  develop 
a  sparse-detector,  high-speed  wavefront  sensor.  To  demonstrate  the  technique,  a 
computational  simulation  of  a  two-dimensional  heated  jet  was  used.  We  explored  the 
feasibility  of  using  multiple,  small-aperture,  probe  beams  placed  along  the  streamwise 
direction  to  use  a  variation  of  Eq.  20  from  Appendix  B,  in  constructing  a  wavefront 
of  a  would-be  large-aperture  collimated  beam  traversing  the  same  flowfield  at  an  in¬ 
stant  in  time.  Our  method  made  use  of  time-series  data  of  the  wavefront  derivatives 
at  the  probe-beam  locations  “measured”  both  prior  to  and  after  the  desired  “con¬ 
struction”  time.  Although  it  appears  to  be  a  straight-forward  exercise,  embedded 
in  the  derivation  of  Eq.  20  is  the  fact  that  the  fluid  structures  that  give  rise  to  the 
“convecting”  optical  aberrations  on  the  wavefront  are  not  frozen  (i.e.,  the  structures 
evolve).6,13  In  addition,  the  convection  velocities  change  with  time  and  from  structure 
to  structure  (more  than  one  of  which’s  aberrating  footprint  is  always  present  in  the 
wavefront).  Fig.  1  (in  the  main  text  of  this  report)  is  a  time  series  of  one-dimensional 
wavefronts  ( X/D  is  the  nondimensionalized  distance  in  the  flow  direction)  due  to 
propagation  through  the  numerical  heated-jet  flowfield  of  Ref.  13.  The  Xj D-vs-time 
slopes  of  the  hills  and  valleys  indicate  the  convection  velocities  of  the  aberrations  on 
the  wavefront,  and  the  amplitude  of  the  OPD  vs.  X/D  gives  the  structure  of  the 
aberrations  on  the  wavefront.  Examination  of  these  features  in  Fig.  1,  demonstrate 
the  facts  that  the  convection  velocities  are  time-varying  and  the  aberration  struc¬ 
tures  evolve  in  time;  thus  the  frozen-flow  hypothesis  inherent  in  Eq.  20  is  not  strictly 
true.  As  described  in  Ref.  13,  cross-correlation  techniques  between  adjacent  probe 
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beams  can  be  used  to  obtain  both  a  running  estimate  of  the  convection  velocities  and 
the  evolution  rates;  these  allow  for  the  construction  of  a  time  series  of  time-resolved 
wavefronts,  and  the  ability  to  place  estimates  on  the  accuracy  of  the  measurements, 
respectively.  Wavefront  sensors  developed  using  these  methods  are  now  referred  to 
as  SABT  (Small-Aperture  Beam  Technique)  wavefront  sensors.  SABT  sensors  have 
demonstrated  a  bandwidth  of  100  kHz17,18  and  have  been  applied  to  optical  wavefront 
measurements  for  propagation  through  heated  jets,6  compressible  shear  layers,18  and 
heated- water  shear  layers.15  To  date  the  SABT  wavefront  sensor  is  the  fastest  wave- 
front  sensor  yet  available;  its  high  bandwidth  is  directly  attributable  to  the  fact  that 
it  is  a  sparse-detector  sensor.  The  sparseness  of  the  sensor  means  that  the  required 
number  of  interrogations  to  acquire  the  information  to  construct  a  full-aperture  wave- 
front  is  reduced  by  one  to  two  orders  of  magnitude  from  that  required  to  obtain  the 
same  resolution  from  a  Shack-Hartmann  wavefront  sensor. 


47 


APPENDIX  D.  SHEAR  LAYER  VELOCITY  FIELD 

SIMULATION 

The  initial,  transitionally-turbulent  flowfield  in  a  shear  layer  is  dominated  by  the 
shear-layer  rollup  caused  by  the  inviscid  Kelvin-Helmholtz  instability  mechanism.41 
This  rollup  has  been  successfully  simulated  using  various  inviscid  and  pseudo-inviscid 
methods.42-45  The  shear  layer  for  the  present  study  was  modeled  using  the  discrete 
vortex  method  (DVM).  DVM’s  are  traditionally  applied  to  incompressible  flows  be¬ 
cause  they  rely  on  the  global  influence  a  discrete  vortex  induces  on  the  flowfield 
through  the  Biot-Savart  law.  In  a  compressible  flow,  the  speed  of  sound,  a,  is  finite 
so  this  influence  would  be  expected  to  travel  along  characteristic  lines.46  Trying  to 
incorporate  compressible  effects  into  the  DVM,  though  perhaps  possible  in  princi¬ 
ple,  would  add  considerable  complexity  and  significantly  increase  the  computational 
demands  of  the  simulation.  On  the  other  hand,  within  certain  Mach- number  re¬ 
strictions,  experimental  observations  of  compressible  shear  layers  (such  as  those  in 
Refs.  34,  37,  and  39)  have  found  similar  vortex  structures  and  mean  velocity  profiles 
to  those  found  in  incompressible  shear  layers.  Thus,  there  must  be  some  validity  to 
applying  DVM  directly  to  a  weakly-compressible  shear-layer  flow.  This  validity  is 
worth  exploring  further. 

D.l.  Evidence  of  Weak  Compressibility 

The  Biot-Savart  Law  is  a  consequence  of  decomposing  the  velocity  V  of  fluid  particles 
in  an  infinitesimal  control  volume  into  rotational  and  irrotational  parts.72  It  can  be 
shown  (see  Refs.  31  and  40,  for  example),  that  the  effect  of  compressibility,  is  to 
introduce  a  dilatation  source  term  A  into  the  irrotational  part  of  the  decomposition. 
This  term  is  ultimately  neglected  in  DVMs  as  they  are  normally  implemented.  Since 
the  mean  velocity  profiles  become  self-similar  for  x  corresponding  to  AEDC  stations  2 
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and  3,  any  mean  portion  of  A  would  vary  only  slowly  in  the  x  direction;  therefore,  we 
would  expect  only  an  optical  "piston”  aberration  or  perhaps  a  slight  tilt  aberration 
which  could  be  neglected.  The  order  of  magnitude  of  the  fluctuating  portion  of 
A  can  be  found  from  the  fluctuating  divergence,  V  •  u',  where  u'  is  the  fluctuating 
portion  of  V  resulting  from  a  Reynolds-type  decomposition.46  According  to  Smits  and 
Dussauge,46  the  ratio  of  the  fluctuating  divergence  to  the  magnitude  of  the  velocity 
gradients  of  the  energy-containing  eddies,  u' / A,  is  given  by 

^«10(7-l)Mt2  (24) 

u  j  A 

where  7  is  the  ratio  of  specific  heats  for  the  gas,  Mt  =  \fu?/a  is  the  turbulent  Mach 
number,  and  overbars  signify  time  averages.  For  the  AEDC  flowfield,  Mt  ~  0.09 
whereby 

«  0.03  «  1  (25) 

v!  j  A 

suggesting  that  the  compressibility  is  weak  and  the  fluctuating  divergence  might  be 
neglected.46,73 

Experimental  condsiderations  also  argue  that  A  has  only  a  second-order  effect 
on  the  AEDC  velocity  field.  Mixing  layers  can  be  characterized  by  the  convective 
Mach  number  Mc  1  of  the  largest  flow  structures.  Mc  1  is  defined  as39 


Mcl  = 


Ui-Uc 

ai 


(26) 


where  U  is  the  velocity  in  the  x-direction,  a  is  the  speed  of  sound,  and  the  convection 
velocity  Uc  is  given  by 

(2?) 

ai  +  a2 

where  1  and  2  refer  to  the  high-  and  low-speed  streams,  respectively.  For  the  con¬ 
ditions  tested  by  Hugo  et  al.  at  AEDC,  the  high-speed  side  of  the  shear  layer  was 
only  Mi  =  0.8;18  this  means  there  were  no  shocks  in  the  flow.  Additionally,  with  a 
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low-speed  side  of  0.1  Mach,  the  convective  Mach  number  Mcl  seen  by  the  largest- 
scale  flow  structures  would  be  only  McX  «  0.33  for  this  case.33-38  Further,  it  has  been 
shown38,37  that  for  compressible  mixing  layers  with  subsonic  Mi ,  the  reduction  in  vi¬ 
sual  spreading  rate  5viz/x  (measured  from  schlieren  photographs)  for  the  compressible 
mixing  layer  relative  to  the  incompressible  case  {8ViZ/x)inc  is  less  than  10%  as  shown 
in  Fig.  20. 
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Figure  20.  Normalized  5^  growth  rate  variation  with  Mc\.  Original  data  compiled 
in  Fig.  9  of  Ref.  37. 


The  shear  layer’s  spreading  rate  can  also  be  defined  in  terms  of  the  growth 
with  increasing  x  of  the  vorticity  thickness  5U  given  by34 


U1-U2 

(dU/dy)m 


(28) 


The  (compressible)  vorticity  thickness  spreading  rate,  8^/x,  for  this  Mcl  actually  re¬ 
mains  unchanged  from  the  incompressible  spreading  rate  as  shown  in  Fig.  21.  More- 
elaborate,  compressible  numerical  simulations  also  have  shown  that  changes  in  shear 
layer  structure  do  not  occur  at  low  Afcl. 24,25,74  Thus,  it  appears  that  the  effect  of  com¬ 
pressibility  might  be  reasonably  neglected  in  the  approximation  of  the  (irrotational 
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part  of  the)  velocity  field  for  the  AEDC  test  conditions. 


Figure  21.  Normalized  5 p  growth  rate  variation  with  Mc i-  Data  as  compiled  in 
Fig.  17  of  Ref.  39. 

D.2.  Velocity  Field  Model 

The  shear  layer  for  the  present  study  was  modeled  using  the  discrete  vortex  method,  a 
technique  first  performed  by  Rosenhead.75  The  method  has  been  subsequently  refined 
using  finite-core  point  (line)  vortices76  and  redistributing  the  vorticity  by  inserting 
additional  vortices  as  the  vorticity-induced  velocity  fields  spread  the  vortex  elements 
apart. 13,77,78  Although  an  inviscid  calculation,  momentum  diffusion  was  modeled  by 
a  temporal  growth  in  the  size  of  the  rotational  cores  of  the  discrete  vortices.22,43  The 
primary  advantage  of  using  a  discrete  vortex  method  is  its  ability  to  model  a  shear 
flow  with  a  fairly  small  computational  domain  thereby  producing  considerable  com¬ 
putational  savings.  Modeling  the  shear  layer  using  multiple  vortex  layers  gives  the 
capability  of  representing  any  arbitrary  velocity  profile  across  the  shear  layer.78  Hugo 
used  only  a  single  vortex  layer  for  each  side  of  his  jet  to  minimize  the  computational 
demands  of  the  model.13  Despite  the  relative  crudeness  of  this  model,  Hugo’s  shear 
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layer  structure  agreed  reasonably  well  with  those  of  Ghoniem’s45  multiple  layer  model. 
Hugo’s  numerical  simulation  produced  reasonable  qualitative  and  quantitative  agree¬ 
ment  with  his  experimental  heated  jet.  A  more  complete  discussion  of  discrete  vortex 
methods  is  given  in  Chapter  4  of  Ref.  22. 

The  velocity  fields  for  all  cases  of  the  present  study  were  produced  by  a  single, 
infinite  shear  layer  modeled  with  discrete  vortices  in  a  manner  similar  to  Hugo.13,22 
The  splitter  plate  and  shear  layer  near  the  plate  trailing  edge  were  modeled  by  discrete 
vortices.  The  influences  of  vortices  far  from  the  plate’s  trailing  edge  were  modeled 
analytically.  The  splitter  plate  vortices’  positions  were  indexed  by  one  temporally- 
consistent  position  which,  in  effect,  kept  the  splitter  plate  configuration  fixed;  once  an 
individual  plate  vortex  was  indexed  past  the  trailing  edge,  it  became  another  shear 
layer  vortex  and  was  allowed  to  convect  with  time  due  to  the  net  induced  velocity  at 
its  position  produced  by  all  the  other  vortices  in  the  model.  The  velocities  on  each 
side  of  the  splitter  plate,  Ui  -  261.04  m/s  and  U2  =  34.7  m/s,  were  chosen  to  match 
those  of  the  AEDC  shear  layer.  These  yielded  an  initial  circulation  density,  jd,  for 
the  vortices  of  jd  =  U2  —  Ui  =  —  AU  =  —226.34  m/s  and  an  imposed  convection 
speed  Ucn  =  {Ui  4-  U2)f 2  =  147.9  m/s.  The  velocity  ratio  across  the  shear  layer  was 
r  =  U2/Ui  =  0.132.  The  initial  vortex  core  diameter  St  for  the  free  vortices  was 
varied  to  model  different  splitter-plate-boundary-layer  thicknesses;  during  this  study, 
5i  =  17.25,  34.5,  and  69  mm.  With  an  initial  core  diameter  and  a  core  growth  rate 
function  based  on  the  growth  rate  of  a  laminar  shear  layer,  the  (instantaneous)  shear 
layers  shown  in  Fig.  22  were  produced. 

The  mean  velocity  profiles  can  be  collapsed  by  suitably  nondimensionalizing 
the  velocity79 

D*-S“5tl  +  eri&r)1  (29) 
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y  position  (m)  y  position  (m)  y  position  (m) 


x  position  (m) 


x  position  (m) 


6i  =  69.00  mm 


5i  —  34.50  mm 


Si  =  17.25  mm 


Figure  22.  DVM  shear  layer  rollup  variation  with  initial  thickness  Si.  Thick  gray 
lines  and  circular  windows  illustrate  vortex  “roller”  size  relative  to  the  AEDC  test 
section.  Smaller,  dotted  boxes  show  fields  used  in  aero-optical  computations. 
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(30) 


and  the  y— position80 

*  V  ~  2/0.5 
y  =  — - - 

where  y0. 5  is  the  y- position  where  U  =  Ucn.  Using  this  nondimensionalization,  the 
computed  mean  velocity  profiles  compare  well  with  both  incompressible  and  weakly 
compressible  (Mc  =  0.51)  experimental  measurements  as  shown  in  Fig.  23.  DVM 
profiles  are  shown  at  several  x-stations;  the  agreement  increases  with  increasing  x. 


Figure  23.  Comparison  of  DVM  mean  velocity  profiles  (<5,/  2  =  34.5  mm)  with  exper¬ 
imental  measurements  in  incompressible  (Ref.  81)  and  weakly-compressible  (Ref.  80) 
shear  layers 

Similarly,  the  rms  streamwise  and  lateral  velocity  fluctuations  as  well  as  the 
Reynolds  stress  can  be  compared  to  the  corresponding  experimental  results  as  shown 
in  Figs.  24  and  25.  In  Figs.  24  and  25,  the  peak  DVM  values  increase  with  increasing 
x.  Similarity  of  the  DVM  velocity  fluctuation  and  Reynolds  stress  profiles  would  be 
expected  to  occur  at  larger  x  than  required  for  mean  velocity  profiles.82  The  variation 
in  profile  magnitude  and  shape  with  x  suggests  that  similarity  conditions  have  not  yet 
been  reached  for  the  <5*  shown.  Similarly-shaped  longitudinal  fluctuation  profiles  were 
obtained  by  Ashurst  in  a  DVM43  and  in  the  water  tunnel  experiments  of  Browand 
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and  his  colleagues.35,83  Winant  and  Browand’s  measurements  showed,  in  fact,  that 
the  magnitude  of  the  central  peak  increased  with  increasing  x  while  the  secondary 
peak  (low-speed  side  “shoulder”)  gradually  disappeared,  taking  on  a  shape  normally 
associated  with  fully-turbulent  shear  layers.35  Such  a  local  maximum  is  predicted  by 
Michalke’s  linear  stability  theory  for  spatially-growing  shear  layers51  and  suggests 
that  the  DVM  behaves  as  if  the  splitter-plate  boundary  layers  were  laminar. 


Figure  24.  Comparison  of  DVM  rms  longitudinal  (left)  and  lateral  (right)  velocity 
fluctuation  profiles  (Si/2  =  34.5  mm)  with  experimental  measurements  in  incompress¬ 
ible  (Ref.  81)  and  weakly-compressible  (Ref.  80)  shear  layers.  A  compressible-Euler 
numerical  result  is  also  shown  (Ref.  84) 

The  growth  of  6, u  with  increasing  x  for  our  numerical  simulation  is  shown  in 
Fig.  26.  The  shear  layer  growth  rate  for  all  initial  core  sizes  is  the  same.  Once  the 
shear  layer’s  thickness,  <5W,  becomes  large  enough  (versus  vortex  core  size),  the  Kelvin- 
Helmholtz  instability  can  be  resolved  and  shear-layer  roll-up  begins.  Upstream  of  this 
point  (for  a  given  initial  vortex  core  diameter),  the  instability  cannot  be  spatially 
resolved  so  5U  ~  Si.  According  to  the  empirical  literature,39  the  shear  layer  growth 
rate  at  the  AEDC  Mc\  can  be  computed  using  Eq.  2  with  Cj  =  0.085  (and  Sw  in 
place  of  Sviz).  For  the  conditions  used  in  our  simulation  (r  =  0.132  and  s  =  1), 
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Figure  25.  Comparison  of  DVM  Reynolds  stress  profiles  (5*/2  =  34.5  mm)  with 
experimental  measurements  in  incompressible  (Ref.  81)  and  weakly-compressible 
(Ref.  80)  shear  layers.  A  compressible-Euler  numerical  result  is  also  shown  (Ref.  84) 

Eq.  28  would  predict  a  growth  rate  of  0.130.  The  growth  rate  for  the  simulated  shear 
layer,  A5w/(x  —  x0)  —  0.139  (with  x0  =  0.06  m)  is  within  7%  of  this  predicted  value. 
As  such,  we  felt  the  numerical  simulation  reasonably  modeled  the  velocity  field  of  a 
weakly-compressible  shear  layer. 


x  position  (m) 

Figure  26.  Growth  of  discrete  vortex  shear  layer  in  terms  of  vorticity  thickness. 
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